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1.0     PURPOSE  OF  THE  STUDY 

Several  considerations  derived  from  the  Virgin  River  adjudication  study  have 
indicated  the  need  of  having  available  a  long-term  series  of  natural  flows  for  the  upper 
reaches  of  the  North  Fork  of  the  Virgin  River  (NFVR).    More  specifically,  at  a  river  site 
located  upstream  from  the  north  boundary  of  Zion  National  Park,  above  the  area  known  as 
Narrows  Canyon. 

To  achieve  this  objective,  a  study  was  implemented,  which  consists  of  the  stochastic 
modeling  and  synthetic  data  generation  of  daily  flows.    Historical  data  of  daily  streamflows 
compiled  by  the  U.S.  Geological  Survey  at  several  locations  in  the  region  were  utilized.   The 
stochastic  modeling  of  daily  flows  was  developed  following  current  modeling  procedures. 

This  report  is  organized  in  the  following  manner.    Section  2.0  defines  the  general 
area  of  interest,  the  study  site  and  specifies  the  hydrological  data  to  be  used.    Section  3.0 
carries  out  a  preliminary  data  analysis,  including  extension  of  flow  records  and  the  estimation 
of  long-term  flow  characteristics  at  the  site  of  interest.    Section  4.0  concentrates  in  the 
analysis  of  the  time  dependence  structure  of  flows  in  the  upper  basin  of  the  NFVR.    Section 
5.0  deals  with  the  identification  of  the  flow  stochastic  model.    And  finally  Section  6.0, 
explains  the  generation  process  of  daily  flows  and  discusses  some  potential  applications  and 
limitations  of  the  synthetically  generated  flows.    To  conclude  this  report,  a  list  of  references 
and  appendices  with  numerical  results  are  included. 

The  reader  will  find  repeatedly  references  to  a  previous  report  by  Diaz  (1992),  which 
provides  background  information  about  the  general  hydrologic  conditions  in  the  region.    The 
report  is  entitled  "Streamflow  Characterization  at  Zion  National  Park,  Utah",  presented  to 
the  National  Park  Service  in  August  1992. 


2.0     STUDY  AREA 

The  general  setting  of  Zion  National  Park  within  the  upper  basin  of  the  Virgin  River 
can  be  found  in  the  previous  report  and  several  other  listed  references.    The  very  brief 
description  provided  in  this  section  concentrates  on  the  North  Fork  of  the  Virgin  River, 
which  runs  through  Zion  National  Park  in  the  North-South  direction.    The  river  area  for  the 
purpose  of  this  study  is  located  in  the  upper  portion  of  the  North  Fork  basin,  outside  the 
Park's  boundaries,  and  upstream  from  the  Narrows  Canyon.    More  specifically,  the  study 
site  is  situated  in  the  proximity  of  a  potential  dam  site  selected  by  the  State  of  Utah,  and 
depicted  in  Figure  1  as  "Bullock  Site"  (approximate  location).  The  cartographic  location  of 
the  Bullock  dam  site  is  identified  by  Township  39S,  Range  9VV  and  Section  32,33. 
According  to  the  Utah  Division  of  Water  Resources,  the  Bullock  site  is  located  in  the  Carmel 
formation  immediately  above  the  Navajo  formation. 

The  NFVR  subbasin  encompasses  a  total  of  123.61  square  miles.    The  stream  extends 
for  approximately  38.8  miles,  from  the  watershed  divide  in  the  Cascades  Spring  area  at 
approximately  elevation  9700  ft,  to  the  lowest  portion  at  the  confluence  with  the  East  Fork  of 
the  Virgin  River  (EFVR),  elevation  3,770  ft.    Water  accretion  from  Cascade  Springs 
contributes  much  of  the  flow  to  the  highest  reaches  of  the  North  Fork.    Deep  snow  packs  and 
summer  thunderstorms  contribute  the  non-uniform  component  of  the  surface  runoff,  whereas 
discharge  from  springs  and  seeps  basically  comprise  the  base  flow. 

The  flow  data  available,  in  the  form  of  mean-daily  values,  has  been  collected  and 
compiled  by  the  U.S.  Geological  Survey.    The  three  USGS  gaging  stations  within  the  area  of 
interest  for  this  study  are  listed  in  Table  1.    The  location  of  these  stations  are  also  shown  in 
Figure  1 ,  referenced  by  their  corresponding  site  number.    Gaging  stations  at  sites  (8)  and  (9) 
are  both  located  in  the  upper  portions  of  the  basin,  in  close  proximity  to  each  other,  only  3.3 
miles  apart,  and  also  close  to  the  study  site.    For  all  practical  purposes,  this  study  assumes 
that  flows  registered  at  site  (9),  USGS  station  No.  09405450  -  NFVR  above  Zion  Narrows 
Canyon,  are  representative  of  natural  flows  conditions  at  the  study  site  (Bullock  site). 


Table  1.    List  of  Gaging  Station  Used  in  this  Study 


Site 

USGS 

Station  Name 

Period  of 

Drainage 

Elevation 

No. 

No. 

and  Location 

Record 

Area  (mr) 

(ft) 

(8) 

09405420 

North  Fork  Virgin  River 
helow  Bullock  Canyon 

1975-1984 

30.0 

6420. 

(9) 

09405450 

North  Fork  Virgin  River 
above  Zion  Narrows 

1979-1984 

42.0 

6000. 

(12) 

09405500 

North  Fork  Virgin  River 
near  Springdale 

1926-1991 

344.0 

3970. 

Fig.l    North  Fork  of  the  Virgin  River  Basin,    Study  Area 


3.0  PRELIMINARY  DATA  ANALYSIS 

3.1  Differences  in  Runoff  Characteristics 

The  stations  listed  in  Table  1  have  observed  mostly  natural  (virgin)  streamflows, 
except  for  some  diversions  upstream  from  the  gage  near  Springdale.   There  are  clear 
differences  between  flow  characteristics  at  the  upper  versus  the  lower  portions  of  the  basin. 
To  illustrate  these  differences,  Figure  2  displays  average  flow  regime  conditions  at  sites  (8) 
and  (12),  for  the  same  10-Yrs  period  of  record,  from  1975  to  1984. 

Graph  2A  shows  average  hydrographs  at  both  locations  comprised  of  mean-daily 
flows  (heavy-solid  lines).  Besides  the  marked  contrast  in  water  yield  as  a  consequence  of  the 
difference  in  the  size  of  the  contributing  drainage  areas,  the  periodic  standard  deviation 
(light-solid  lines)  also  shows  marked  differences.    Flow  variability  is  also  shown  in  Graph  2B 
but  expressed  in  the  form  of  the  coefficient  of  variation  (Cv),  computed  as  the  ratio  of  the 
standard  deviation  of  daily  flows  divided  by  their  mean  value.    Notice  the  more  uniform  (less 
variance)  flow  pattern  at  site  (8)  in  comparison  to  that  at  site  (12).    The  more  uniform  flow 
pattern  at  high  elevations  is  attributed  to  the  predominant  relatively  uniform  types  of  water 
sources  contributing  runoff,  such  as  snowmelt  producing  runoff  and  groundwater  contribution 
to  surface  flow.    In  contrast,  the  site  near  Springdale  located  at  low  elevation  in  the  NFVR 
basin,  displays  a  flow  pattern  where  besides  the  two  water  sources  mentioned  above,  a  third 
source,  rainfall,  plays  a  very  important  role. 

The  surface  runoff  pattern  is  a  result  of  the  precipitation  regime  in  the  region. 
Precipitation  frequency  distribution  at  Zion  national  Park  (Station  #  9717)  displays  a  bimodal 
shape,  with  the  two  peaks  during  the  months  of  February  and  August,  and  with  the  lowest 
frequency  during  the  month  of  June.    The  first  peak  corresponds  to  the  winter  snowfall  and 
the  second  one  to  summer  thunderstorms.    Sudden  produce  sharp  flow  increases  of  relatively 
short  duration.    They  can  take  place  all  year  around,  but  mostly  outside  the  snowmelt  season. 
Thunderstorms  are  particularly  evident  during  the  late  spring  and  summer  seasons  at  the 
lower  elevations. 


600 


390 


Days  (Water  Year,  Oct  1  -  Sep  30) 


Fig. 2    Runoff  Characteristics  at  High  and  Low  Areas  of  the  NFVR  Basin,  [1975-1984]  - 
(A)   Periodic  Mean  and  Standard  Deviation,   (B)  Periodic  Coefficient  of  Variation 


3.2       Extension  of  Flow  Record 

To  improve  the  reliability  of  the  periodic  mean  and  variance  of  daily  flows  estimated 
at  the  study  site  (9),  its  historical  period  of  record  is  extended  based  on  data  at  neighboring 
stations,  sites  (8)  and  (12). 

As  indicated  in  Table  1,  10  water-years  of  daily  flows  were  collected  at  site  (8),  from 
1975  to  1984.   The  station  downstream,  site  (9),  has  only  6  years  of  daily  flows  available, 
from  1979  to  1984.   The  close  proximity  of  the  two  stations  and  the  six  years  of 
simultaneous  daily  records  give  us  the  opportunity  to  extend  the  record  of  flows  at  site  (9), 
from  six  to  ten  years,  in  a  reliable  manner.   The  shortest  record  is  extended  based  on  the 
values  registered  at  the  station  with  the  longest  sequence  during  the  non-concurrent  period  of 
record.    The  procedure  is  as  follows: 

1.         First,  daily  flows  were  aggregated  and  regressed  seasonally,  i.e.,  corresponding  to  12 
monthly  seasons.    For  instance,  the  186  (31  days  x  6  years)  daily  flow  values 
corresponding  to  the  month  of  October  at  site  (9)  are  regressed  versus  the 
corresponding  186  values  at  site  (8).   Flow  values  were  initially  transformed  to  the 
logarithmic  domain  in  order  to  more  closely  resemble  a  Gaussian  process.    Second 
order  polynomials  were  selected  as  the  most  suitable  relationship  between  the  two 
variables,  outperforming  the  linear  model  in  all  12  seasons.   The  adopted  second 
order  regression  model  takes  the  form, 

Y  =  a+bX+cX2   +    ^(1  ~r2)  oye  EQ-C1) 

where:  Y  :  dependent  variable 

X  :  independent  variable 

a,b,c  :  parameters  of  the  model 

r  :  cross-correlation  between  series  X  and  Y 

<ry  :  standard  deviation  of  series  Y 

e  :  random  normal  variable  with  mean  '0'  and  variance  '1' 


Table  2.1  shows  the  basic  statistics  of  daily  flows  (after  the  Log-transformation  and 
monthly  aggregation)and  the  cross  correlation  value  of  the  corresponding  series, 
which  is  a  dimensionless  measure  of  the  correlative  association  between  the  two 
series.   The  right  portion  of  Table  2. 1  contains  the  parameters  of  the  regression 
model  and  the  coefficient  of  determination  r2,  which  indicates  the  proportion  of  the 
variance  of  Y  explained  by  the  regression  function.   The  twelve  sets  of  regressed 
values  and  their  corresponding  fitted  models  were  plotted  as  shown  in  Appendix  I. 
Notice  that  in  order  to  explain  all  the  variance  of  the  dependent  variable,  a  random 
component  has  been  added  to  the  regression  term.    In  this  manner,  the  extra  variance 
(1-r2),  not  explained  by  the  pure  regression  term,  is  introduced  in  the  right-hand  side 
of  Eq.(l). 

2.         Flows  at  sites  (8)  and  (9)  were  also  correlated  at  the  monthly  and  yearly  levels, 

which,  as  expected,  show  an  even  higher  degree  of  association  than  at  the  daily  level. 
Results  of  the  linear  models  adopted  at  the  monthly  and  yearly  levels  are  shown  in 
Tables  2.2  and  2.3  respectively.    Estimated  monthly  and  yearly  flows  at  site  (9) 
during  the  period  1975-1979,  based  on  flows  measured  at  site  (8)  during  that  same 
period,  were  used  to  adjust  the  regressed  daily  values  at  site  (9),  such  that  monthly 
volumes  first,  and  annual  volumes  second,  were  also  preserved. 

In  summary,  as  a  result  of  the  record  extension  procedure  described  above,  a  total  of 
ten  water-years  of  mean-daily  flows  became  available  at  site  (9),  NFVR  above  Narrows 
Canyon,  from  1975  to  1984. 
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4.0  LONG-TERM  FLOW  CHARACTERISTICS  AT  THE  STUDY  SITE 

4.1  Estimation  of  Mean  Annual  Discharge 

In  general,  it  is  accepted  that  given  a  period  of  flow  record  extending  for  at  least  10 
years,  the  analyst  can  compute  reasonably  good  estimates  of  the  long-term  flow 
characteristics  at  the  site.    Unfortunately,  this  should  not  be  considered  the  case  at  our  study 
site.    Reviewing  Figure  6  from  the  previous  report,  which  shows  annual  discharge  for  the 
NFVR  near  Springdale  for  the  period  of  record  [1926-1991],  it  can  be  inferred  that  the  10- 
yrs  subperiod  from  1975  to  1984,  as  an  average,  it  would  indicate  a  water  availability 
considerably  larger  than  what  it  should  be  expected  in  the  long  run.   The  same  information  is 
repeated  below  in  Figure  3,  where  the  subperiod  [1975-1984]  has  been  separated  from  the 
remaining  record  for  visualization  purposes. 
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Fig. 3   Annual  Discharge  at  the  NFVR  Near  Springdale,  Site  (12) 
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Notice  that  during  those  ten  years,  the  two  wettest  (1980  and  1983)  and  the  driest 
(1977)  water- years  on  record  took  place.   The  average  annual  discharge  for  the  10-Yrs, 
95,010  Ac-Ft,  is  a  27.5%  larger  than  the  65-Yrs  average,  74,500  Ac-Ft,  at  site  (12).   Then 
it  can  be  hypothesized  that  similar  flow  conditions  existed  in  the  upper  portion  of  the  basin, 
and  consequently,  that  the  mean  annual  discharge  computed  from  the  short  record  at  station 
09405450  is  not  representative  of  the  long-term  flow  conditions  at  site  (9). 

To  confirm  the  hypothesis  above,  annual  discharge  for  the  NFVR  at  sites  (9)  and  (12) 
were  correlated  for  the  10-Yrs  period  with  concurrent  data  (recall  earlier  extension  of  flow 
record  at  site  (9)  from  6  to  10  Yrs.).   The  10  pairs  of  data  and  the  adopted  linear  regression 
function  are  shown  in  Figure  4.    As  expected,  annual  flows  measured  near  Springdale  and 
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Fig. 4   Regression  of  Annual  Flows  at  Site  (9)  versus  Site  (12),  Period  [1975-1984], 
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above  Narrows  Canyon  show  a  high  linear  correlation,  yielding  a   coefficient  of 
determination  r2  of  0.95.   This  high  value  of  r2  is  even  more  meaningful  when  considering 
that  the  correlated  values  cover  the  whole  range  of  annual  flows  ever  registered  at 
Springdale,  from  the  driest  to  the  wettest  water-years.    The  linear  relationship  in  Figure  4 
can  therefore  be  utilized  to  estimate  the  long-term  mean  annual  discharge  for  the  NFVR 
above  Narrows  Canyon,  computed  as  10,860  Ac-Ft  (15.03  cfs),  based  on  the  measured  65- 
Yrs  mean  annual  discharge  near  Springdale  of  74,500  Ac-Ft. 

4.2       Estimation  of  Periodic  Mean  and  Periodic  Standard  Deviation 

Once  the  long-term  mean  annual  discharge  at  site  (9)  was  estimated,  the  following 
step  was  to  estimate  the  long-term  periodic  mean  and  periodic  standard  deviation  of  daily 
flows.   The  procedure  adopted  for  this  task  was  similar  to  the  one  adopted  to  transfer 
information  from  site  (8)  to  site  (9).    But  this  time,  the  correlation  of  daily  flows  was  carried 
out  between  sites  (9)  and  (12).    Site  (12)  is  the  closest  location  to  site  (8)  with  a  relatively 
long  period  of  record  of  daily  flows,  62  years. 

In  the  same  way  done  earlier  in  Section  3.2,  after  estimating  a  complete  month  of 
daily  values  (by  means  of  the  daily  regression  equations),  daily  flows  were  adjusted  (when 
necessary)  to  match  the  estimated  flow  value  at  the  monthly  level.    Similarly,  once  a  whole 
year  of  daily  flows  were  generated,  they  were  also  subject  to  adjustment  in  order  to  conform 
with  the  estimated  discharge  at  the  annual  level.    As  a  result  of  this  sequential  flow 
adjustment  procedure,  the  extended  database  of  daily  flows  at  site  (9),  from  1926  to  1988 
(1927  excluded),  yields  an  average  mean  annual  discharge  practically  equal  to  the  value 
found  in  Figure  4  ,  10,860.  Ac-Ft  (or  its  equivalent  15.26  cfs). 

To  help  visualize  the  change  in  average  annual  discharge  as  a  result  of  the  length  of 
the  period  of  analysis  being  considered,  Figure  5  displays  three  different  average  annual 
hydrographs.    Each  one  of  them  represents  average  daily  flows  based  on  the  three  subperiods 
discussed  so  far,  [1979-1984],  [1975-1984]  and  [1927-1988],  the  last  one  also  named  "long- 
term".    The  hydrographs  shown  in  Figure  5  are  the  smooth  representation  of  the  actual 
periodic  functions,  after  being  fitted  by  Fourier  series. 
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Fig. 5   Average  Annual  Hydrographs  at  Site  (9)  as  a  Function  of  the  Period  of  Analysis 


Notice  that  as  the  length  of  the  period  of  analysis  increases,  the  mean  annual 
discharge,  i.e.,  the  area  under  the  curve,  decreases.    This  is  consistent  with  the  findings  from 
Figure  3,  which  shows  extremely  wet  water-years  during  the  6-Yrs  subperiod,  still  wet 
conditions  for  the  10-Yrs  subperiod,  and  the  approximation  to  long-term  (normal)  conditions 
when  the  period  of  62-Yrs  is  being  considered.    Moreover,  the  hydrographs  become  smother 
when  going  from  the  shorter  to  the  longer  period  of  record  as  a  result  of  less  sampling 
variability. 

Another  important  aspect  when  analyzing  the  structure  of  daily  runoff  is  the 
variability  of  daily  flows  around  its  average  value.  To  compare  the  degree  of  flow 
variability  between  measured  and  regressed  values,  the  periodic  coefficient  of  variation  (Cv) 
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was  computed  from  both  series.    There  were  not  substantial  differences  between  the  Cv 
values  computed  from  the  10-Yrs  and  62-Yrs  subperiods,  although  the  longer  series  shows 
slightly  higher  values  during  the  summer  season.    These  higher  values  can  be  attributed  to 
the  effect  of  summer  thunderstorms  registered  at  site  (12),  that  to  some  extent,  are  carried 
along  by  the  flow  regression  procedure  into  site  (9). 

In  summary,  either  one  of  the  periodic  Cv  functions  provides  reasonably  good 
estimates  of  the  variance  of  flows,  and  therefore,  can  be  used  to  estimate  the  long-term 
standard  deviation  at  Site  (9).    Figure  6  shows  the  adopted  long-term  sample  mean,  sample 
standard  deviation  and  their  corresponding  smoothed  periodic  functions  for  the  NPVR  site 
above  Narrows  Canyon.    The  numeric  values  of  the  periodic  functions  shown  in  Figure  6 
were  included  as  Appendix  II. 
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5.0  DEPENDENCE  STRUCTURE  OF  DAILY  FLOWS 

The  understanding  of  the  dependence  structure  of  short  interval  time  series  is  crucial 
for  the  stochastic  modeling  of  daily  time  series.    The  periodicity  in  parameters  that  describe 
or  measure  the  dependence  of  successive  values  of  daily  flows,  k-lag  apart,  can  be  studied 
from  a  physical  and  a  statistical  point  of  view.    The  latter  is  the  main  subject  of  investigation 
in  this  section. 

The  time  dependence  structure  of  daily  runoff  is  a  highly  site-specific  flow 
characteristic.    For  that  reason,  the  10-Yrs  period  of  record  [1975-1984]  at  Site  (9)  was 
adopted  to  study  the  dependence  structure  of  daily-flows  in  the  upper  region  of  the  NFVR 
basin.    The  10-Yrs  period,  although  short,  it  should  provide  a  good  representation  of  flow 
sequencing  due  to  the  alternation  of  very  dry,  very  wet  and  normal  water-years  during  that 
particular  period. 

5.1  Normalization  and  Standardization  of  Daily  Flows 

In  studying  the  properties  of  autocorrelative  dependence  of  flows  decisions  should  be 
made  on  whether  to  study  the  structure  of  the  original  flow  series,  or  alternatively, 
"standardized"  and/or  "normalized"  transformations  of  the  flow  series.    For  instance,  when 
analyzing  the  periodic  skewness  and  kurtosis  coefficients  of  daily  flows  at  site  (12),  see 
Figures  13  and  14  in  the  previous  report,  it  becomes  apparent  that  the  marginal  distributions 
of  daily  flows  change  from  one  day  to  another,  yielding  a  (periodic)  highly  non-stationary 
process.    Moreover,  the  values  of  the  skewness  and  kurtosis  coefficients  depart  substantially 
from  zero  and  three  respectively,  indicating  how  different  these  distributions  are  from  the 
normal  distribution.   This  typical  departure  from  normality  of  daily  flows  is  also  observed  at 
site  (9),  making  the  structural  models  available,  developed  for  Gaussian  processes, 
unapplicable  for  the  case  of  daily  flows. 

A  way  of  circumventing  the  problem  of  lack  of  normality  of  daily  flows  is  to  find  an 
appropriate  transformation  of  the  asymmetric  series,  such  that  the  transformed  series  is  closer 
to  normality.   Typically  in  hydrology,  the  logarithmic,  the  power  and  the  Wilson-Hilferty 
functions  are  used  to  normalize  daily  flows.   The  normalized  series  Yv  T  is  obtained  as 
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*"».,  ■  /x<*,»>  Eq-(2) 

where  the  function  /,  is  chosen  based  on  the  skewness  test  of  normality. 

On  the  other  hand,  standardization  of  a  series  removes  the  within  the  year  seasonal 
pattern  of  its  means  and  standard  deviations.    The  standardized  series  is  computed  by 

=     (^-ffl,)  ^  (3) 

sx 

for  y=l,2,...N  with  N  the  number  of  years,  and  r=l,2,...,w  with  w  =  365  the  number  of 
time  intervals  per  year.    Equation(3)  yields  the  standardized  variable  Zu>  T  where  m  T  and  s  T 
are  the  periodic  means  and  standard  deviation  respectively  of  the  series  being  standardized. 
This  transformation  restores  stationarity  in  the  two  parameters,  yielding  a  transformed 
process  with  approximately  mean  zero  and  variance  one.    Third  and  fourth  order  moments 
remain  unaffected  after  the  standardization.    Depending  on  the  characteristics  of  the  flow 
series  at  hand,  standardization  is  applied  before  or  after  normalization.    Table  3  presents 
results  from  all  the  attempted  transformations  of  daily  flows  at  site  (9)  during  the  period 
[1975-1984].    Results  are  presented  as  average  for  the  365  values  of  the  periodic  mean, 
periodic  standard  deviation,  periodic  skewness  and  periodic  kurtosis  coefficients  of  the 
resulting  series  after  the  transformation  .    First  case  implies  no  transformation  (original 
domain). 

Table  3   Average  Statistics  After  Transforming  Daily  Flows  Series  at  Site  (9) 

Type  of  Average  value  of  periodic  ...  #  Days  Rejected 

Transformation  m,  sr  g,  kr  99%  C.L. 


Y„ 

r    =    X  „,  , 

19.34 

14.92 

0.94 

4.61 

148 

Y, 

t  =  Login  (X„,  r) 

1.12 

0.36 

-0.55 

4.76 

128 

Y. 

>r  =  Log_e  (X„,r) 

2.58 

0.83 

-0.55 

4.76 

128 

Y, 

>r  =  Sqrt  (X„,r) 

4.00 

1.54 

0.30 

4.03 

37 

Y„ 

,  T  =  Log  Wilson-H 

ilferty  (X  „, ,) 

0.00 

1.01 

0.01 

4.11 

35 
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Table  3  indicates  the  Log  Wilson-Hilferty  (LWH)  as  the  most  suitable  transformation, 
yielding  the  smallest  number  of  days  rejecting  the  normality  hypothesis,  with  35days  outside 
the  99%  confidence  limits.    Unlike  the  other  transformations  where  the  series  is  first 
normalized  and  then  standardized,  LWH   involves  first  transforming  the  original  series  into 
the  logarithmic  domain  as  in  Eq.(2).    Then,  the  series  is  standardized  using  Eq.(3),  where 
mr  and  sr  are  computed  in  the  Log-domain.    Finally,  the  Wilson-Hilferty  transformation  is 
carried  out  using  Eqs.(4)  and  (5)  to  arrive  at  Zu  T .    Matalas  (1967)  gives 


U  ,  T 


G,(y) 


G,(y)  y\. 


+  1 


1/3 


-  1 


ct(y) 


Eq(4) 


valid  for  values  of  the  skewness  Gr(y)^0.    If  GT(y)=0  no  transformation  is  necessary.    The 
Fourier  fitted  series  of  the  daily-skewness  coefficients  in  the  log(e)  domain  is  used  for  the 
LWH  transformation.    Furthermore,  Y*„  T  is  given  by  McGinnis  and  Sammons  (1970) 


Y* 


O,  T 


max[7b>t/-2/GT(y)] 
min  [7U(t,  -2  /Gx(y)] 


if  Gx{y)  >0 
if  G,(y)  <0 


Eq.(5) 


After  completing  the  LWH  transformation,  normality  and  stationarity  of  the 
standardized  series  is  reached  for  all  practical  purposes.    Figure  7(A)  shows  the  periodic  first 
two  moments  of  the  LWH  transformed  series.    Notice  that  the  periodic  patterns  were  totally 
removed  from  both  series,  with  mean  zero  (average  )  and  variance  one  (average) 
respectively,  as  a  result  of  the  standardization  .    Also,  the  LWH  transformation  produced 
average  coefficients  of  skewness,  Figure  7(B),  and  kurtosis  which  are  close  to  zero  and  three 
respectively,  as  expected  for  the  normal  marginal  distributions,  regardless  of  their  small 
variations  around  the  averages.    Due  to  the  relatively  small  number  of  years  of  record 
available  for  the  analysis  and  the  large  number  of  time  intervals  under  consideration,  the 
Fourier  series  fitting  of  the  periodic  characteristics  of  the  historical  series  is  desirable.    The 
parametric  method  of  normalization  allow  us  to  decrease  the  number  of  parameters  of  the 
periodic  components  and  to  smooth  out  sampling  variability. 
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Fig. 7    Periodic  Statistics  of  Daily-Flows  at  Site  (9)  in  the  Log  Wilson-Hilferty  Domain. 
(A)  Mean  and  Standard  Deviation,   (B)  Skewness  Coefficients 
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5.2       Analysis  of  Dependence  Structure 

The  seasonal  serial  correlation  coefficients  r  k  r  of  daily  flows  at  site  (9)  is  displayed 
by  Figure  8.   The  periodic  function  rk  T  is  estimated  by  the  expression 


k,x 


N 

I 

u=l 


2  (*».*  ~m^  (*•.*♦*  ~mx*k) 


11/2 


E  (X».t-m<)2   £  (*„.*♦*-*,♦*> 


u*l 


u-1 


Eq.(6) 


where  all  variables  have  been  previously  defined.   The  most  important  case  is  for  k=l,  the 
periodic  first  autocorrelation  coefficient.   The  r^  T  function  definitely  shows  much  less  within 
the  year  variation  than  the  equivalent  periodic  function  computed  for  the  NFVR  near 
Springdale,  as  shown  in  Figure  15  of  the  previous  report. 
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Fig. 8    Lag-1  Autocorrelation  Function  of  Standardized  Daily  Flows  at  Site  (9) 

The  small  variance  shown  by  r,  T  is  displayed  mostly  outside  from  the  part  of  the 
year  in  which  the  snowmelt  is  the  major  runoff  contributor  (r=  180-280).  Notice  that  the 
majority  of  the  values  are  contained  within  a  relatively  narrow  range  (0.8-0.9). 
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Furthermore,  the  periodicity  in  r,  T  can  be  considered  negligible  in  practical  terms.    In  other 
words,  it  can  be  safely  assumed  that  the  lag-1  serial  autocorrelation  of  daily  flows  at  site  (9) 
do  not  significantly  vary  from  day-to-day,  and  therefore,  it  can  be  represented  by  a  uniform 
value  of  r,  r  ■  0.86   for  r  =  1  ,2,  ...,  365.   Besides  r,  r,  autocorrelations  functions  for 
other  lags  were  also  computed  (though  not  shown  in  Figure  8).    rk  T  functions  for  k=  1,2  and 
3  resulted  almost  undistinguishable  among  themselves.   This  finding  give  us  an  early 
indication  of  the  uniformity  and  persistence  of  mean-daily  flows  at  high  elevation  watersheds 
of  the  NFVR  hydrologic  region. 

Because  of  the  negligible  periodicity  in  the  serial  autocorrelation  of  flows,  the 
seasonal  time  series  Zu >T  can  be  lumped  into  a  non-periodic  series  ZT,  with  r=l,2  ...  N', 
where  N'=N  .  w.    For  the  case  under  consideration,  N'  =  10  Yrs  x  365  days/Yrs  =  3650 
occurrences  of  the  Z  variable.    Consequently,  the  k*  autocorrelation  coefficient  is  computed 
as  a  function  of  lag  (k)  only 

rk  =  f(k)  Eq.(7) 

The  computation  of  rk  can  still  be  carried  out  by  Eq.(6),  used  to  compute  rk  T ,  except 
that  the  time  index   r   is  dropped.    Figure  9  shows  the  rk  function,  known  as  the 
correlogram,  for  the  standardized  series  of  daily  flows  at  site  (9)  in  the  original  flow-domain, 
which  uses  the  mean  and  standard  deviation  of  the  entire  Zv  r  series.    Note  the  very  smooth 
and  monotonic  decay  of  the  correlogram  from  r,  to  r390,  with  values  greater  than  zero  even  at 
the  end  of  the  year.   The  slow  decay  rate  indicates  a  hydrologic  system  with  very  long  water 
storage  memory.    The  positive  time  dependence  of  rk  also  explains  the  large  persistence  of 
flows  of  the  same  type,  that  is,  low  values  mostly  followed  by  low  values,  and  conversely, 
high  values  followed  by  high  values.   The  correlogram  shows  no  indication  of  any  remaining 
periodicity  in  the  standardized  series.    The  r,  value  of  the  correlogram  is  0.956.   The  same  r, 
value  computed  in  the  Log  Wilson-Hilferty  domain  drops  slightly  to  0.938. 

To  contrast  the  difference  in  flow  dependence  structure  between  the  upper  and  the 
lower  regions  of  the  NFVR  basin,  the  correlogram  of  daily  flows  at  site  (12),  during  the 
same  10- Yrs  period,  was  also  included  in  Figure  9.    Notice  the  much  faster  decay  of  rk  at 
site  (12),  particularly  during  the  first  10  lag-periods.   The  larger  influence  of  precipitation  in 
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Fig. 9   Correlogram  of  Standardized  Daily  Runoff  at  Sites  (9)  and   (12)   -  Flow  Domain 

the  form  of  rainfall  at  the  lower  elevations  is  the  main  factor  responsible  for  the  difference  in 
the  correlation  structure  of  the  residual  series  between  the  two  areas. 

Another  tool  for  analyzing  flow  dependence  structure  is  the  partial  autocorrelation 
function  (PAF),  also  displayed  in  Figure  9,  for  site  (9).   PAF  measures  the  degree  of  linear 
association  between  successive  autocorrelation  coefficients.   Together  with  the  correlogram, 
help  the  analyst  in  the  identification  of  the  order  of  the  linear  dependence  model  that  more 
closely  will  resemble  the  flow  structure  of  the  residual  series.   The  first  eight  PAF 
coefficients  are  significant  at  the  95  percent  probability  levels,  given  approximately  by 
+2A/N'  =  0.033.   Then,  it  oscillates  within  the  confidence  band.   These  results  would 
suggest  a  stochastic  model  with  a  relatively  large  number  of  autorregresive  terms. 
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5.3       Identification  of  the  Flow  Generation  Model 

As  indicated  in  Section  5.1,  the  periodicity  in  the  mean  and  standard  deviation  of 
daily  flows  was  removed  by  using  a  Fourier  series  representation  of  YT  and  sT  and  by 
creating  the  standardized  series  Zu  T .    By  removing  the  periodicity  in  the  mean  and  the 
standard  deviation,  the  series  Zv  T  became  second-order  stationary,  provided  the  serial 
autocorrelation  function  can  be  considered  approximately  stationary  as  demonstrated  in 
Figure  8. 

The  typical  flow  pattern  of  the  NFVR,  with  low  flows  mainly  from  groundwater 
contribution  and  higher  flows  produced  by  rainfall  events  and  snowmelt,  demands  a  flow 
generation  model  that  it  can  accommodate  this  mixed  behavior.    Models  with  autoregressive 
(AR)  and  moving  average  (MA)  components,  known  as  ARMA  models,  are  recommended 
under  this  conditions.   More  specifically,  ARMA  models  with  constant  parameters  are 
expected  to  perform  satisfactory  modeling  mean-daily  flows,  provided  that  the  flow 
correlation  function  is  approximately  stationary  (Salas  et  al.,  1985).   The  mixed 
autoregressive-moving  average  model  for  fitting  the  periodic  hydrologic  series  Y„  T  is 
expressed  by 

yut  =  mx  *  sxZViX  Eq.(8) 

where  Zu  r  will  be  represented  by  an  ARMA(p,q)  model  with  (p)  autoregressive  parameters 
and  (q)  moving  average  parameters,  constant  in  time, 


Z,  =  f  *,  zt_,  -  f  8,  e^  +  eT  Eq.(9) 

J-i  7^1 

where  <f>it  j  =  l,2,...p  and  6it  i=l,2,...q  are  the  coefficients  of  the  model  to  be  evaluated 
from  data,  and  er  is  an  independent  normal  variable.   The  visual  display  of  the  original  series 
and  the  behavior  of  the  autocorrelation  function  coupled  with  that  of  the  PAC  have  revealed 
that  a  relatively  large  number  of  <£'s  and   0's  parameters  will  be  required. 

The  estimation  of  the  ARMA(p,q)  parameters  is  a  complex  process.    Details  of  the 
estimation  procedure  are  not  provided  here,  however  they  can  be  found  in  Salas  et  al.(1985, 
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pp.223).    Preliminary  estimates  of  the  parameters  of  tentative  models  were  obtained  by  the 
Method  of  Moments,  and  a  posteriori,  more  accurate  estimates  of  the  parameters  were 
obtained  by  the  Non-linear  Least-square  algorithm.    Computations  were  performed  by  a 
computer  program.    Analysis  of  a  parameter  estimation  output  is  analyzed  in  Section  5.4. 

Several  competing  ARMA(p,q)  models  were  fitted  to  the  standardized  ZT  series.    The 
Akaike  Information  Criterion  (AIC),  Akaike  (1974),  a  test  used  to  select  the  best  among 
competing  models  was  utilized  to  select  the  orders  (p)  and  (q).   The  AIC  is  computed  by 

AIC  =  N  Ln(oe2)  +  2  (p+qr)  Eq.(lO) 

where  N  is  the  sample  size  and  o;  is  the  variance  of  the  estimated  residuals  series  obtained 
after  fitting  the  prospect  ARMA(p,q)  model  to  the  data  series.    The  model  that  has  the 
minimum  AIC  value  is  the  theoretical  choice.    The  model  identification  procedure  started 
with  low  order  models,  for  instance  an  ARM A(  1,1),  and  it  continued  progressively 
overfitting  the  time  series  until  the  minimum  AIC  value  was  obtained.    At  the  same  time,  the 
goodness-of-fit  of  the  prospect  model  was  verified.   The  analysis  showed  values  of  AIC  (and 
o})  relatively  similar  for  models  with  AR  terms  greater  or  equal  than  three.    However,  only 
an  ARMA(4,2)  satisfied  the  the  Parsimony  Conditions.    Consequently,  an  ARMA(4,2)  model 
was  adopted  for  further  analysis. 

5.4       Diagnostic  Checking 

After  estimating  the  parameters  of  the  identified  model,  it  is  necessary  to  verify  the 
validity  of  the  underlying  assumptions.    The  following  program  output  provides  the  values  of 
the  parameters  4>j  and  B{ ,  checks  the  parsimony  of  the  model  parameters  and  tests  the 
residuals  series  for  independence  and  normality.  All  these  requirements  are  derived  from  the 
estimation  procedures  that  are  commonly  used  in  structural  analysis  of  hydrologic  processes. 
The  details  of  these  tests  may  be  found  in  Salas  et  al.(1985).    All  conditions  were  accepted 
except  the  hypothesis  of  normality  of  the  residuals  eT ,  rejected  by  the  skewness  test  of 
normality.    Daily  streamflows  typically  display  this  departure  from  Gaussian  processes. 
However,  in  this  particular  case,  the  results  of  the  normality  test  can  be  ignored  since  the 
normality  requirement  of  the  transformed  series  was  already  verified  (after  LWH,  Table  3). 
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Table  4.    Parameter  Estimation  and  Diagnostic  Check  for  ARMA(4,2)  Model 


Project 

Name  of  River  Basin 

USGS  gaging  Station 

Type  of  Data 

Data  Units 

Period  of  Analysis 


Zion  National  Park 

North  Fork  of  Virgin  River 

09405450  -  NFVR  above  Narrows  Canyon 

Mean  Daily  Flows 

No  Units,  Variable  Standardized  and  Normalized 

using  Log_Wilson-Hilferty  Transformation 

1975-1984 


Statistical  Characteristics  of  the  Dependent  Variable 


>  Type  of  series 

>  Mean  of  the  Series 

>  Standard  Deviation 

>  Skewness  Coefficient 


Standardi zed/Normal i zed (LWH) /Stationary 
-0.001 
0.979 
0.252 


Sample  Autocorrelation  Function 
k    -0.95CL    r(k)   +0.95CL 


1 

-0.032 

0.938 

0.032 

2 

-0.032 

0.912 

0.032 

3 

-0.032 

0.899 

0.032 

4 

-0.032 

0.884 

0.032 

5 

-0.032 

0.873 

0.032 

6 

-0.032 

0.867 

0.032 

7 

-0.032 

0.862 

0.032 

k 

-0.95CL 

r(k) 

+0.95CL 

8 

-0.032 

0.856 

0.032 

9 

-0.032 

0.849 

0.032 

10 

-0.032 

0.840 

0.032 

50 

-0.032 

0.734 

0.032 

100 

-0.032 

0.659 

0.032 

200 

-0.032 

0.536 

0.032 

300 

-0.032 

0.346 

0.032 

Parameter  Estimation  of  AR  and  MA  Parameters 

>  Moment  Estimators  >  Least  Squares  Estimators 

Par.  No.    AR     MA  Par.  No.    AR      MA 


1 

1.000 

0.387 

2 

0.000 

0.128 

3 

0.000 

4 

0.000 

1 

1.091 

0.505 

2 

0.157 

0.350 

3 

-0.147 

4 

-0.103 

Residual  Variance: 0. 103 


Residual  Variance:0.099 


Stationarity  and  Invertibility  Check  for  ARMA  Models  (Parsimony  Test) 
>  AR  Parameters  >  MA  Parameters 

Par. No.      AR  Value        Root  Par. No.    AR  Value        Root 


1.096 
0.188 
■0.207 
•0.080 


2.234 
1.979 
1.913 
1.145 


Stationarity  Condition  is  Accepted 


0.501 
0.361 


-1.709 
-1.673 


..Invertibility  Condition  is  Accepted 
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A  Portmanteau  Lack  of  Fit  Test  for  Independence  of  Residuals 


Computed  Critical  Value 
Chi-square  Tabulated  Value 
Degrees  of  Freedom 
Confidence  Level 


834.9 
976.9 
906 
0.950 


Hypothesis  of  Independence  is  Accepted 


Skewness  Test  for  Normality  of  Residuals 


Tabulated  Test  Value 
Computed  Test  Value 
Probability  Level 


0.067 
2.180 
0.950 


Hypothesis  of  Normality  is  Rejected 
Recommendation  ==>  Transform  Series 


Overfitting  Check  Using  the  Akaike  Information  Criterion 


AIC 

Of 

ARMA 

(4,2) 

Model 

-8428. 

7 

AIC 

Of 

ARMA 

(5,2) 

Model 

-8428. 

9 

AIC 

Of 

ARMA 

(4,3) 

Model 

-8405. 

7 

AIC 

Of 

ARMA 

(3,2) 

Model 

:   -8398. 

7 

AIC 

Of 

ARMA 

(4,1) 

Model 

-8421. 

6 
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6.0  GENERATION  OF  SYNTHETIC  FLOWS 

6.1  Data  Generation  Scheme 

The  flow  generation  process  is  started  by  using    Eq.(9)  recursively  to  generate 
synthetic  values  of  ZT.    The  computational  process  is  initiated  by  arbitrarily  assigning    (p) 
initial  values  of   'Z'    and  computing    (q)    initial  values  of   V  .    The  first  generated  50 
values  of   ZT  are  discarded  to  avoid  the  transient  effect  of  the  initial  conditions.    Rewriting 
Eq.(9)  according  to  the  identified  ARMA(4,2)  model,  the  following  expression  results 

z<  -  2>;zw  "  EM-i  +  /0T099  $T  Eq.(ll) 

i=i  1=1 

with  parameters  <j>-t  ,  0,  and  a(  obtained  from  Table  4.    £T  is  the  standardized  normal 
variable  with  mean  zero  and  variance  one,  generated  using  the  Box-Muller  formula  written  as 

e1  =    [2  Ln{l/ur)  ]  1/2  cos  (2ttu2)  Eq  (12) 

e2  =    [2  Ln(l/ux)  }  in  sin(27ru2) 

where  u,  and  u2  are  two  independent  uniformly  distributed  (0,1)  random  numbers.  A 
subprogram  that  generates  normal  random  numbers  was  made  part  of  the  flow  generation 
model. 

Having  obtained  the  sequence  of  ZUT,  the  backward  VVilson-Hilferty  transformation  is 

applied  using 


Y.  _2_   K(y) 

U'T  '      Gr(y) 


GAy) 


Gx{y) 


Eq.(13) 


in  which  Gr(y)?^0,  where   Gr(y)   is  the  skewness  coefficient  of  the  daily  series  in  the 
Log(e)  domain.    The  variable    Y*UiT  is  further  transformed  by 


26 


max[Y\z,  -2/GT(y)  ]     ...    ifGT(y)>0 

min[ir*U(t,-2/Gt(y)]    ...    ifGx(y)<0 

rUjT  =  zU/T  ...   if  GT(y)  =0 


r«>.t  -  \  min[y*U(t,-2/Gt(y)]    ...   ifGt(y)<0  Eq.(14) 


Yu<  r  becomes  the  generated  data  in  the  Log-domain  (theoretically  with  mean  zero  and 
variance  one).    Finally,  the  data  Xu  T  (flows  in  the  original  domain)  is  obtained  from 

XViX  =  exp[FT  +  st(y)  yv#t]  Eq.(15) 

where  Y,.and    sT(y)  are  the  estimated  long-term  periodic  mean  and  standard  deviation  of  the 
logarithms  of  daily  flows  for  the  NFVR  above  Narrows  Canyon.    Equivalent  periodic 
functions  for  the  extended  period  [1926-1988],  but  in  the  original  flow  domain,  were  shown 
in  Figure  6.    A  computer  program  was  implemented  for  synthetic  data  generation. 

6.2      Validation  of  Generated  Data 

The  last  step  in  any  stochastic  modeling  effort  relates  to  the  validity  of  the  chosen 
model.    In  order  to  evaluate  the  reliability  of  the  fitted  ARMA(4,2)  to  model  daily-flows  at 
site  (9)  some  statistical  analysis  will  be  performed  based  on  actual  data  generation  as 
described  in  Section  6.1.    Before  concentrating  in  the  analysis,  it  should  be  pointed  out  that 
stochastic  models  are  built  to  "resemble"  the  main  statistical  characteristics  of  the  historic 
series  in  the  statistical  sense.    That  is,  it  should  not  be  expected  that  synthetic  series  based  on 
the  model  will  give  "exactly"  the  same  statistical  characteristics  as  depicted  by  the  historic 
series,  but  rather,  we  should  look  for  a  close  resemblance  to  those  characteristics. 

Fifteen  hydrological  series  of  daily  flows,  70  years  long  each  (for  a  total  of  1050 
years),  were  generated.    Basic  statistics  were  computed  and  compiled  separately  for  each 
hydrological  series  in  the  original  flow  domain.    Plotted  in  Figure  10   are  the  long-term  basic 
statistics,  daily-means  in  Graph  (A)  and  daily-standard  deviation  in  Graph  (B)  (darker  lines), 
together  with  the  upper  and  lower  confidence  bands  of  the  generated  statistics.    The 
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10   Historical  Daily  Means  -  Graph  (A),  and  Standard  Deviation 
Confidence  Bands  of  Generated  Statistics 


Graph  (B),  and 
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confidence  band  is  computed  as  positive  and  negative  one-standard  error  relative  to  the 
average  of  the  15  series. 

In  general,  the  results  show  that  long-term  statistics  are  satisfactorily  reproduced. 
Notice  that  the  generated  statistics  follow  very  closely  to  the  periodic  pattern  of  the  historic 
series.    The  confidence  bands  of  the  generated  statistics  enclose  the  historical  statistics, 
although  in  the  case  of  the  standard  deviation,  the  latter  lies  close  to  the  lower  bound  during 
the  snowmelt  season.    This  is  not  unusual  since  the  band  is  only  one  standard  error  above 
and  below  the  generated  mean.    Presumably,    the  larger  than  expected  flow  variability  of  the 
synthetically  generated  series  is  due  to  the  increase  in  flow  variance  introduced  at  site  (9) 
during  the  flow  extension  procedure,  see  Section  4.2. 

The  generated  first  serial  autocorrelation    function  r, >T  (not  shown)  resulted  slightly 
smoother  than  the  historical  one,  since  a  constant  (non-periodic)  autocorrelation  coefficient 
was  adopted  during  the  generation.    Moreover,  the  average  value  of  the  r,  T  function  (for  r  = 
1  to  365)  increased  from  0.86  for  the  historical  series  (see  Figure  8)  to  0.92  for  the  synthetic 
series.    The  flow  dependence  structure  of  the  synthetic  series  was  also  compared  with  the 
historic  series  by  plotting  the  correlograms.    Figure  1 1  displays  the  two  corresponding 
correlograms,  showing  a  very  good  reproduction  of  the  correlation  structure  of  flows. 
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Fig.  1 1    Periodic  Correlogram  of  Historical  and  Synthetic  Series  of  Daily  Flows 
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This  flow  dependence  aspect  of  the  whole  modeling  effort  is  considered  critical  when 
judging  the  performance  of  the  stochastic  model.    Particularly  in  this  case,  due  to  the 
extremely  large  persistence  of  flows    in  the  upper  region  of  the  NFVR  basin.    Note  that  the 
correlogram  shown  in  Figure  1 1  at  site  (9)  is  the  same  as  the  one  shown  in  Figure  9  (darker 
line),  except  that  in  the  case  of  Figure  9  the  flow  series  was  standardized  prior  to  applying 
Eq.(6),  and  consequently,  removing  the  periodicity  embedded  in  the  daily  flow  series. 

In  summary,  it  can  be  concluded  that  the  resemblance  of  generated  and  historical 
statistics  is  as  expected  and  desired,  making  the  model  useful  for  the   generation  of  synthetic 
flows  in  the  area  of  the  NFVR  above  Zion  Narrows.    However,  it  is  important  to  point  out 
that  the  stochastic  model  developed  in  this  study  has  the  characteristic  of  being  univariate, 
that  is,  the  stochastic  component  ot  the  model  was  developed  using  information  from  a  single 
river  site.    While  this  is  not  a  limitation  for  using  the  model  as  a  tool  to  simulate  surface 
runoff  conditions  at  a  local  level,  it  prevents  the  analyst  from  using  the  generated  flows  in 
conjunction  with  other  existing  flow  records  in  the  area,  in  a  more  regional  type  of  analysis. 
In  case  that  several  concurrent  (in  time)  hydrologic  series  were  necessary  for  the  analysis, 
the  multivariate  stochastic  analysis  and  multivariate  modeling  will  become  imperative. 
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Appendix  I. 


Regression  of  Mean-Daily  Flows  at  Site  (9)  versus  Mean-Daily  Flows 
at  Site  (8) 
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Appendix  II  : 


Long-Term  Estimates  of  Periodic  Mean  and  Periodic  Standard 
Deviation    of  Dailv  Flows  for  the  NFVR  above  Narrows  Canvon 


Day 

. .  E  s  t  i  ma  t  ed 

...  Fitted 

Day 

. .  E  s  t  i  ma  t  ed 

...  Fitted 

No. 

Mean 

Std.Dev 

Mean   Stc 

I.Dev 

No. 

Mean 

Std.Dev 

Mean   Std 

.Dev 

1 

11.72 

10.05 

10.13 

6.43 

61 

11.45 

6.52 

12.73 

6.47 

2 

12.65 

10.21 

10.24 

6.51 

62 

13.29 

6.05 

12.81 

6.40 

3 

10.20 

5.62 

10.35 

6.58 

63 

12.70 

5.51 

12.88 

6.33 

4 

10.60 

5.85 

10.46 

6.64 

64 

13.30 

6.42 

12.95 

6.26 

5 

10.32 

5.43 

10.58 

6.69 

65 

12.04 

6.11 

13.01 

6.19 

6 

10.96 

4.20 

10.69 

6.73 

66 

12.17 

5.05 

13.07 

6.11 

7 

11.14 

5.18 

10.80 

6.76 

67 

11.29 

5.90 

13.11 

6.04 

8 

9.25 

4.46 

10.91 

6.78 

68 

12.23 

5.69 

13.15 

5.97 

9 

10.69 

5.90 

11.01 

6.78 

69 

13.94 

6.28 

13.18 

5.91 

10 

10.03 

5.47 

11.11 

6.77 

70 

13.66 

6.94 

13.20 

5.85 

11 

11.05 

5.49 

11.20 

6.74 

71 

14.74 

7.24 

13.21 

5.80 

12 

9.27 

6.61 

11.29 

6.70 

72 

12.57 

5.20 

13.21 

5.77 

13 

11.15 

6.39 

11.37 

6.65 

73 

13.19 

6.20 

13.20 

5.74 

H 

11.55 

6.33 

11.45 

6.59 

74 

13.15 

6.16 

13.18 

5.72 

15 

10.90 

6.18 

11.52 

6.51 

75 

12.48 

5.87 

13.14 

5.72 

16 

11.66 

6.55 

11.58 

6.43 

76 

13.30 

6.02 

13.09 

5.73 

17 

12.17 

7.07 

11.63 

6.34 

77 

12.67 

5.58 

13.03 

5.75 

18 

10.99 

6.20 

11.68 

6.24 

78 

11.71 

5.71 

12.96 

5.78 

19 

12.63 

6.63 

11.72 

6.14 

79 

13.44 

6.81 

12.88 

5.83 

20 

12.23 

6.59 

11.75 

6.03 

80 

12.96 

6.53 

12.79 

5.88 

21 

11.04 

6.75 

11.77 

5.93 

81 

12.57 

6.13 

12.69 

5.95 

22 

10.59 

6.22 

11.79 

5.83 

82 

12.30 

5.48 

12.58 

6.02 

23 

11.35 

5.58 

11.80 

5.73 

83 

14.11 

6.85 

12.46 

6.11 

24 

11  .01 

6.15 

11.80 

5.63 

84 

13.43 

6.39 

12.33 

6.19 

25 

11.09 

5.58 

11.80 

5.54 

85 

12.13 

5.98 

12.19 

6.29 

26 

12.82 

6.75 

11.80 

5.47 

86 

13.71 

6.66 

12.05 

6.38 

27 

12.23 

5.11 

11.79 

5.40 

87 

12.74 

6.70 

11.90 

6.48 

28 

11.52 

4.90 

11.77 

5.34 

88 

12.81 

6.72 

11.75 

6.57 

29 

11.89 

5.19 

11.75 

5.29 

89 

11.68 

6.70 

11.60 

6.66 

30 

12.21 

5.90 

11.73 

5.26 

90 

12.98 

7.72 

11.45 

6.75 

31 

11.35 

5.43 

11.71 

5.24 

91 

11.98 

6.17 

11.29 

6.83 

32 

12.05 

5.08 

11.69 

5.23 

92 

12.50 

6.52 

11.14 

6.91 

33 

11.21 

6.48 

11.67 

5.24 

93 

8.62 

5.59 

10.98 

6.97 

34 

11.75 

5.36 

11.64 

5.27 

94 

8.61 

5.76 

10.84 

7.02 

35 

12.72 

5.66 

11.62 

5.30 

95 

9.85 

6.51 

10.69 

7.07 

36 

11.08 

4.76 

11.61 

5.35 

96 

8.57 

5.95 

10.56 

7.10 

37 

11.77 

5.03 

11.59 

5.41 

97 

8.99 

6.13 

10.43 

7.12 

38 

11.68 

5.34 

11.58 

5.48 

98 

9.98 

6.26 

10.30 

7.13 

39 

12.28 

5.18 

11.57 

5.56 

99 

9.33 

5.82 

10.19 

7.13 

40 

11.87 

5.07 

11.57 

5.64 

100 

10.37 

6.68 

10.09 

7.12 

41 

11.47 

4.68 

11.57 

5.74 

101 

9.39 

6.76 

10.00 

7.10 

42 

11.65 

7.11 

11.58 

5.83 

102 

9.94 

6.31 

9.92 

7.07 

43 

11.49 

5.89 

11.59 

5.93 

103 

9.94 

5.68 

9.85 

7.03 

44 

12.32 

6.72 

11.61 

6.03 

104 

9.82 

6.79 

9.80 

6.99 

45 

13.01 

6.65 

11.64 

6.13 

105 

10.95 

7.72 

9.76 

6.95 

46 

12.00 

5.56 

11.67 

6.22 

106 

10.14 

10.77 

9.74 

6.91 

47 

12.87 

6.45 

11.71 

6.31 

107 

10.05 

7.91 

9.73 

6.87 

48 

12.84 

6.26 

11.76 

6.39 

108 

9.17 

6.65 

9.73 

6.84 

49 

12.31 

6.18 

11.81 

6.46 

109 

9.11 

6.44 

9.75 

6.81 

50 

11.96 

6.69 

11.87 

6.53 

110 

10.84 

10.01 

9.78 

6.80 

51 

11.59 

5.96 

11.93 

6.58 

111 

11.61 

9.05 

9.83 

6.79 

52 

12.30 

6.79 

12.00 

6.62 

112 

11.21 

7.78 

9.89 

6.80 

53 

11.03 

6.25 

12.08 

6.65 

113 

10.39 

7.96 

9.96 

6.83 

54 

12.06 

6.68 

12.15 

6.67 

114 

10.18 

7.13 

10.05 

6.88 

55 

12.53 

7.09 

12.23 

6.67 

115 

10.19 

6.60 

10.14 

6.96 

56 

12.37 

6.70 

12.31 

6.67 

116 

10.37 

7.06 

10.25 

7.05 

57 

11.79 

6.29 

12.40 

6.65 

117 

10.51 

6.89 

10.37 

7.17 

58 

11.97 

6.40 

12.48 

6.62 

118 

11.39 

6.78 

10.49 

7.31 

59 

11.18 

6.32 

12.57 

6.58 

119 

12.93 

7.97 

10.63 

7.49 

60 

13.18 

6.42 

12.65 

6.53 

120 

10.76 

6.93 

10.77 

7.68 

35 


(cont inuat  ion) 


Day 

. .    E  s  t  i  ma  t  ed 

...    Fitted   ... 

Day 

. .    Estimated 

...    F  i  1 1  ed   

No. 

Mean 

Std.Dev 

Mean 

Std.Dev 

No. 

Mean 

Std.Dev 

Mean 

Std.Dev 

121 

12.60 

12.16 

10.91 

7.90 

181 

24.18 

18.39 

23.87 

14.84 

122 

10.86 

9.51 

11.06 

8.15 

182 

25.59 

17.37 

24.49 

15.64 

123 

11.38 

9.32 

11.22 

8.42 

183 

20.92 

12.77 

25.11 

16.47 

124 

10.77 

8.72 

11.37 

8.71 

184 

22.14 

13.59 

25.73 

17.32 

125 

10.90 

8.04 

11.53 

9.03 

185 

19.53 

11.93 

26.35 

18.18 

126 

10.52 

7.37 

11.69 

9.35 

186 

24.98 

14.97 

26.97 

19.05 

127 

10.79 

7.62 

11.84 

9.70 

187 

22.77 

16.68 

27.58 

19.92 

128 

10.90 

7.49 

11.99 

10.05 

188 

24.35 

18.93 

28.18 

20.79 

129 

11.00 

7.19 

12.14 

10.41 

189 

30.44 

24.09 

28.77 

21.64 

130 

10.09 

6.23 

12.29 

10.77 

190 

30.16 

27.24 

29.35 

22.48 

131 

11.47 

7.59 

12.42 

11.13 

191 

25.87 

21.91 

29.91 

23.30 

132 

11.85 

6.21 

12.56 

11.48 

192 

27.39 

20.58 

30.45 

24.09 

133 

10.87 

5.42 

12.68 

11.83 

193 

28.21 

18.86 

30.97 

24.84 

134 

12.15 

6.50 

12.80 

12.16 

194 

31.83 

20.02 

31.46 

25.57 

135 

12.62 

6.99 

12.91 

12.47 

195 

35.53 

22.88 

31.93 

26.25 

136 

12.57 

6.79 

13.02 

12.76 

196 

27.99 

20.86 

32.37 

26.89 

137 

14.11 

27.22 

13.12 

13.02 

197 

28.12 

22.68 

32.78 

27.49 

138 

14.74 

25.27 

13.20 

13.26 

198 

35.40 

31.73 

33.15 

28.03 

139 

12.91 

17.24 

13.29 

13.46 

199 

37.83 

35.59 

33.49 

28.54 

140 

12.86 

18.95 

13.36 

13.62 

200 

32.52 

27.91 

33.80 

28.99 

141 

13.57 

21.72 

13.43 

13.75 

201 

38.06 

30.99 

34.07 

29.39 

142 

14.33 

21.33 

13.49 

13.83 

202 

36.71 

29.30 

34.30 

29.74 

143 

13.93 

26.00 

13.55 

13.87 

203 

40.68 

39.66 

34.49 

30.04 

144 

13.73 

16.44 

13.60 

13.87 

204 

38.83 

37.24 

34.65 

30.29 

145 

13.09 

13.44 

13.65 

13.83 

205 

43.52 

39.83 

34.76 

30.49 

146 

14.98 

11.87 

13.70 

13.75 

206 

45.38 

33.57 

34.83 

30.65 

147 

13.66 

10.65 

13.75 

13.62 

207 

36.72 

26.70 

34.86 

30.77 

148 

15.00 

10.18 

13.80 

13.46 

208 

33.73 

26.68 

34.85 

30.84 

149 

13.22 

8.03 

13.85 

13.26 

209 

42.61 

39.52 

34.80 

30.88 

150 

11.67 

7.96 

13.90 

13.02 

210 

34.57 

31.05 

34.72 

30.88 

151 

12.68 

8.12 

13.96 

12.76 

211 

45.79 

39.99 

34.59 

30.84 

152 

15.97 

9.12 

14.03 

12.47 

212 

33.24 

27.96 

34.42 

30.78 

153 

15.06 

8.24 

14.10 

12.16 

213 

26.95 

22.40 

34.22 

30.68 

154 

14.22 

10.34 

14.18 

11.83 

214 

27.07 

21.87 

33.98 

30.56 

155 

15.74 

9.54 

14.28 

11.49 

215 

27.05 

22.74 

33.71 

30.42 

156 

13.92 

8.10 

14.39 

11.15 

216 

27.16 

24.18 

33.41 

30.25 

157 

16.28 

12.45 

14.51 

10.81 

217 

28.77 

31.87 

33.07 

30.07 

158 

13.20 

8.97 

14.65 

10.48 

218 

30.50 

32.25 

32.71 

29.87 

159 

15.84 

8.89 

14.81 

10.16 

219 

30.38 

29.85 

32.32 

29.65 

160 

15.16 

8.21 

14.99 

9.85 

220 

32.86 

33.48 

31.90 

29.42 

161 

15.89 

9.60 

15.19 

9.58 

221 

28.30 

29.64 

31.47 

29.18 

162 

14.86 

8.69 

15.40 

9.33 

222 

32.62 

34.28 

31.01 

28.93 

163 

18.83 

11.99 

15.65 

9.12 

223 

27.91 

26.86 

30.53 

28.66 

164 

15.81 

9.81 

15.91 

8.96 

224 

27.91 

26.55 

30.03 

28.38 

165 

15.95 

11.51 

16.20 

8.84 

225 

31.44 

30.02 

29.52 

28.10 

166 

17.86 

10.75 

16.51 

8.77 

226 

27.58 

26.22 

29.00 

27.80 

167 

16.72 

9.35 

16.85 

8.76 

227 

28.04 

25.88 

28.47 

27.49 

168 

17.71 

9.43 

17.21 

8.80 

228 

27.58 

26.77 

27.93 

27.18 

169 

17.29 

8.97 

17.60 

8.91 

229 

26.94 

26.10 

27.39 

26.85 

170 

16.63 

8.11 

18.01 

9.08 

230 

28.55 

27.67 

26.84 

26.51 

171 

20.13 

11.21 

18.44 

9.31 

231 

27.05 

27.31 

26.29 

26.15 

172 

16.72 

8.62 

18.90 

9.60 

232 

25.93 

24.47 

25.74 

25.78 

173 

20.94 

12.38 

19.38 

9.96 

233 

22.67 

24.28 

25.19 

25.40 

174 

20.26 

11.22 

19.88 

10.38 

234 

26.64 

28.28 

24.64 

25.01 

175 

21  .41 

11.57 

20.41 

10.86 

235 

23.05 

23.28 

24.10 

24.60 

176 

21.87 

12.87 

20.95 

11.40 

236 

21.79 

21.72 

23.57 

24.17 

177 

20.57 

12.97 

21.51 

11.99 

237 

20.98 

21.36 

23.04 

23.73 

178 

23.12 

15.75 

22.08 

12.64 

238 

20.41 

21.41 

22.52 

23.27 

179 

23.40 

15.81 

22.67 

13.33 

239 

21.37 

21.61 

22.02 

22.79 

180 

27.38 

20.55 

23.27 

14.07 

240 

18.89 

20.39 

21.52 

22.30 
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(cont  inuat  ion) 


Day 

. .  Estimated  . . 

...  Fitted  ... 

Day 

..  Estimated  .. 

...  Fitted 

No. 

Mean 

Std.Dev 

Mean 

Std.Dev 

No. 

Mean 

Std.Dev 

Mean   Stc 

l.Dev 

241 

18.20 

18.93 

21.03 

21.80 

301 

11.25 

6.67 

10.99 

7.65 

242 

18.13 

19.53 

20.56 

21.28 

302 

12.41 

6.99 

11.05 

7.69 

243 

18.21 

18.31 

20.10 

20.74 

303 

9.77 

5.85 

11.11 

7.74 

244 

24.72 

25.29 

19.65 

20.20 

304 

10.01 

6.04 

11.17 

7.80 

245 

26.62 

26.32 

19.22 

19.64 

305 

11.82 

7.52 

11.22 

7.85 

246 

20.99 

21.58 

18.80 

19.07 

306 

13.24 

8.37 

11.28 

7.92 

247 

21.73 

22.82 

18.39 

18.50 

307 

14.10 

8.99 

11.34 

7.99 

248 

21.14 

20.45 

18.00 

17.92 

308 

12.58 

7.64 

11.39 

8.06 

249 

19.94 

19.39 

17.62 

17.34 

309 

12.67 

8.21 

11.44 

8.13 

250 

20.02 

18.90 

17.25 

16.76 

310 

11.06 

10.43 

11.49 

8.21 

251 

20.57 

18.71 

16.90 

16.18 

311 

10.42 

7.00 

11.53 

8.29 

252 

20.11 

17.50 

16.56 

15.60 

312 

12.07 

9.04 

11.57 

8.37 

253 

18.90 

14.80 

16.24 

15.03 

313 

11.06 

7.62 

11.60 

8.44 

254 

18.03 

13.87 

15.92 

14.47 

314 

10.73 

10.16 

11.62 

8.52 

255 

16.62 

13.44 

15.62 

13.92 

315 

12.28 

8.10 

11.63 

8.59 

256 

14.96 

12.51 

15.33 

13.39 

316 

16.17 

12.44 

11.64 

8.66 

257 

14.22 

11.30 

15.06 

12.87 

317 

11.43 

9.25 

11.64 

8.72 

258 

13.44 

10.13 

14.79 

12.37 

318 

9.37 

6.45 

11.63 

8.77 

259 

14.78 

11.10 

14.53 

11.88 

319 

10.04 

10.07 

11.62 

8.82 

260 

13.06 

9.99 

14.29 

11.42 

320 

10.51 

8.72 

11.59 

8.85 

261 

11.18 

8.45 

14.05 

10.99 

321 

12.45 

17.20 

11.56 

8.87 

262 

11.33 

8.12 

13.82 

10.58 

322 

11.69 

16.28 

11.52 

8.88 

263 

13.36 

9.61 

13.60 

10.19 

323 

10.70 

7.38 

11.47 

8.88 

264 

9.74 

6.82 

13.39 

9.83 

324 

9.40 

5.87 

11.41 

8.86 

265 

11.47 

8.17 

13.19 

9.50 

325 

12.63 

8.25 

11.34 

8.83 

266 

13.46 

9.40 

13.00 

9.19 

326 

9.61 

6.30 

11.27 

8.79 

267 

10.67 

7.86 

12.81 

8.92 

327 

10.62 

6.99 

11.19 

8.72 

268 

11.17 

7.81 

12.63 

8.66 

328 

12.35 

8.98 

11.10 

8.65 

269 

12.20 

8.27 

12.46 

8.44 

329 

10.71 

9.02 

11.01 

8.56 

270 

10.39 

7.80 

12.30 

8.24 

330 

10.63 

7.63 

10.91 

8.45 

271 

9.25 

6.67 

12.14 

8.06 

331 

9.74 

7.38 

10.81 

8.33 

272 

9.56 

7.07 

11.99 

7.91 

332 

10.28 

7.08 

10.70 

8.20 

273 

9.42 

6.70 

11.85 

7.78 

333 

10.10 

7.01 

10.59 

8.06 

274 

11.37 

8.61 

11.72 

7.67 

334 

9.35 

6.45 

10.48 

7.91 

275 

10.98 

8.52 

11.59 

7.57 

335 

9.73 

7.69 

10.37 

7.75 

276 

12.93 

8.73 

11.47 

7.50 

336 

10.75 

6.97 

10.26 

7.58 

277 

12.92 

9.06 

11.36 

7.44 

337 

10.13 

6.75 

10.15 

7.41 

278 

9.90 

6.70 

11.25 

7.39 

338 

8.82 

5.41 

10.04 

7.24 

279 

11.53 

7.60 

11.16 

7.36 

339 

9.05 

5.29 

9.94 

7.07 

280 

12.25 

8.71 

11.07 

7.33 

340 

9.14 

5.49 

9.84 

6.90 

281 

10.27 

7.90 

10.99 

7.32 

341 

10.64 

6.03 

9.74 

6.73 

282 

13.72 

10.20 

10.91 

7.31 

342 

10.92 

7.23 

9.65 

6.57 

283 

11.96 

9.44 

10.85 

7.31 

343 

8.96 

7.16 

9.57 

6.41 

284 

12.67 

10.33 

10.79 

7.31 

344 

13.34 

9.77 

9.49 

6.27 

285 

13.17 

10.58 

10.74 

7.32 

345 

8.73 

8.74 

9.42 

6.13 

286 

13.01 

10.72 

10.71 

7.33 

346 

8.43 

5.08 

9.36 

6.01 

287 

10.31 

8.79 

10.67 

7.34 

347 

10.32 

5.73 

9.31 

5.90 

288 

8.92 

6.95 

10.65 

7.35 

348 

10.57 

6.77 

9.26 

5.81 

289 

13.23 

8.07 

10.64 

7.36 

349 

8.60 

5.54 

9.23 

5.73 

290 

12.07 

7.09 

10.63 

7.38 

350 

9.42 

5.38 

9.21 

5.67 

291 

10.56 

5.53 

10.63 

7.39 

351 

8.00 

4.76 

9.20 

5.63 

292 

10.21 

5.37 

10.64 

7.41 

352 

10.08 

5.38 

9.20 

5.60 

293 

10.11 

5.65 

10.66 

7.43 

353 

11.40 

6.05 

9.20 

5.59 

294 

9.93 

7.02 

10.68 

7.44 

354 

9.42 

5.36 

9.22 

5.59 

295 

9.41 

7.40 

10.71 

7.46 

355 

10.23 

5.52 

9.25 

5.61 

296 

9.16 

5.69 

10.74 

7.49 

356 

9.04 

5.22 

9.29 

5.65 

297 

9.88 

6.01 

10.79 

7.51 

357 

9.13 

5.50 

9.34 

5.69 

298 

10.12 

5.77 

10.83 

7.54 

358 

9.67 

6.28 

9.40 

5.75 

299 

11.81 

6.62 

10.88 

7.57 

359 

10.41 

6.90 

9.47 

5.82 

300 

9.60 

5.30 

10.93 

7.61 

360 

8.42 

4.46 

9.55 

5.90 
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(continuation) 


Day 

..  Estimated  .. 

...  Fitted 

NO. 

Mean 

Std.Dev 

Mean 

Std 

.Oev 

361 

8.43 

5.11 

9.63 

5.98 

362 

9.42 

9.21 

9.72 

6.07 

363 

10.51 

6.74 

9.81 

6.16 

364 

10.26 

6.88 

9.91 

6.25 

365 

8.50 

6.91 

10.02 

6.34 

Ave. -> 

15.26 

14.45 

15.26 

14.45 
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